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We consider the capture of inelastic dark matter in white dwarves by inelastic spin-independent 
scattering on nuclei. We show that if the dark matter annihilates to standard-model particles then, 
under the assumption of primordial globular cluster formation, the observation of cold white dwarves 
in the globular cluster M4 appears inconsistent with explanations of the observed DAM A/LIBRA 
annual modulation signal based on spin-independent inelastic dark matter scattering. Alternatively 
if the inelastic dark matter scenario were to be confirmed and it was found to annihilate to standard- 
model particles then this would imply a much lower dark matter density in the core of M4 than 
would be expected if it were to have formed in a dark matter halo. Finally we argue that cold 
white dwarves constitute a unique dark matter probe, complementary to other direct and indirect 
detection searches. 



I. INTRODUCTION 

Many models of dark matter (DM), in particular those 
models where the DM abundance is a consequence of 
thermal freeze-out, require the annihilation of DM into 
standard-model particles which usually implies a non- 
zero DM-nucleon interaction cross-section. These two 
properties can be constrained by a variety of observa- 
tions. The annihilation cross-section and products can 
be constrained by the requirement that DM particles are 
not overproduced during thermal freeze-out and/or by 
limits on fluxes of standard model particles that would 
arise from DM annihilation in the galactic halo. The 
DM-nucleon cross-section can be constrained by direct 
detection experiments which search for the recoil energy 
deposited when DM scatter off nuclei. Limits on a com- 
bination of the DM-nucleon cross-section and the partic- 
ular annihilation products can be placed by considering 
the flux of neutrinos that would result from DM capture 
and subsequent annihilation in the Sun or Earth. 

It has also been pointed out in [T] |2] that as white 
dwarves (hereafter WDs) have no internal energy source, 
and many cold WDs have now been observed, it is pos- 
sible to set limits on the DM-nucleon cross-section by 
showing that the energy released by the annihilation of 
DM particles in their core could contribute significantly 
to the luminosity of the star. 

An interesting possibility for DM-nucleon interactions 
is that the DM could predominantly scatter off nucleons 
inelastically to an excited energy state, so-called 'inelas- 
tic dark matter' (iDM). This idea has been proposed as 
a possible explanation |2{5] for the annual modulation 
signal observed by the DAMA collaboration [S| |7| . The 
key feature in iDM is the mass splitting, S, between DM 
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particles that scatter off nuclei. For iDM explanations of 
the DAMA/LIBRA annual modulation 6 is of the order 
30 — 130 keV (depending on DM mass and couplings). 

The main physical consequence of the inelastic split- 
ting is that the minimum velocity for a DM particle to 
scatter off a nucleus and impart an energy Eji is in- 
creased. This can severely weaken the sensitivity of di- 
rect detection experiments as the number of particles in 
the halo with a large enough velocity to scatter can be 
very small, or even zero in some cases. As a consequence 
the allowed DM-nucleon cross-sections can be orders of 
magnitude larger than for elastic scattering. 

Limits on iDM from capture in the Sun [8l-[T0] are 
promising as the escape velocity of the Sun, Vesc ~ 
600 — 1300 km s~^, is large enough to provide sufficient 
energy to in-falling DM particles to overcome inelastic 
splittings 5 ~ 100 keV. Although these limits are ham- 
pered by the difficulty of detecting neutrinos, progress in 
neutrino telescope exposure means strong limits exist on 
particular models of iDM, for example, limits on sneu- 
trino iDM have been studied in [llj. However, models 
of iDM that annihilate predominantly to e~^e~, fj,'^fj,~, 
77, light hadrons or gluons are immune to these limits as 
these particles either stop before decaying and produc- 
ing neutrinos or don't produce neutrinos at all. In these 
cases the energy deposited in the Sun from DM annihi- 
lations is swamped by the internal energy due to fusion. 
Therefore models of the class described in [I^ can evade 
limits from solar capture. 

The converse is true for limits on capture in WDs as 
it is the deposited energy that is used to set the lim- 
its, and annihilation to neutrinos will only deposit a 
small amount of energy in the WD. Therefore models of 
DM that attempt to explain the DAMA annual modula- 
tion observation with the iDM mechanism, in particular 
models that also aim to offer an explanation of recent 
PAMELA [13^ and Fermi LAT |T4j observations through 
DM annihilation, are subject to limits from capture in 
WDs. 
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In this paper we investigate limits on the iDM sce- 
nario by considering the temperature and luminosity of 
recently observed WDs in our closest globular cluster M4. 
Throughout we assume that the DM annihilates to stan- 
dard model particles. 



II. WHITE DWARVES IN M4 

White dwarves are compact objects made up of a de- 
generate electron core comprised almost entirely of car- 
bon and oxygen. The electron degeneracy prevents any 
contraction and the temperature of this core is too low 
to ignite nuclear fusion reactions. As a result, WDs have 
no internal energy source and release only the thermal 
energy of the non-degenerate ions in the core. For these 
reasons the evolution of WDs can be described as a cool- 
ing process and the age of a globular cluster such as M4 
can be estimated by observing a cut-off at low magni- 
tudes in the WD cooling sequence. This has motivated 
observations of WDs in globular clusters down to very 
low magnitudes. 

Recently this low magnitude cut-off has been observed 
in the globular cluster M4 [T3] . We use the best-measured 
set of data from these observations, subject to all of the 
selection processes detailed in [TS], including the use of 
proper motion measurements, to produce a decontami- 
nated sample. 

We take the data in the form of magnitudes in the 
F606W and F775W Hubble Space Telescope filters and 
convert the colour meoeiv — to an effective tem- 

perature under the assumption that the WD is radiat- 
ing as a black-body. To do this we numerically shine a 
black-body spectrum through the HST filter transmission 
curves, convert magnitudes to the ACS/WFC Vega-mag 
system |16) and correct for reddening and extinction as 
detailed in Section 9 of [H]. We then use the mgoew 
magnitude, the M4 distance modulus, and the effective 
temperature of each star and compare with Vega [53] to 
calculate the luminosity for each star. This data is plot- 
ted in Fig. ^ 

Having obtained the luminosity and temperature it is 
possible to calculate the radius of each WD (under the 
assumption of black-body radiation). Then, using the 
Salpeter equation of state [T7], we calculate an approxi- 
mate mass-radius relationship, which can be used to de- 
termine an approximate mass for each star. In Fig. [5] we 
show the stars in the mass-luminosity plane. 



III. DARK MATTER IN M4 

In order to set limits on the iDM-nucleon cross-section 
it is necessary to estimate the DM density surrounding 
the WD's at the center of M4. Although it is currently 
impossible to do this to with any great accuracy, recent 
developments in the observation and simulated evolution 
of globular clusters embedded in galactic halos now allow 



an estimate of the DM content. Some time ago Peebles 
[l8] suggested that globular clusters may be formed in 
subhalos of DM before falling into galactic halos. How- 
ever observations of order ~ 1 mass-to-light ratios [T5] 
and the tidal stripping of stars from some globular clus- 
ters [54] suggest a significant DM component cannot re- 
side within or without the observed stellar distribution 
[20| . These observations set an upper limit on the DM 
content of globular clusters. 

Recent simulations have shed light on how these results 
can be reconciled with a primordial scenario of globular 
cluster formation through the process of tidal stripping. 
In fact, the presence of globular clusters has been sug- 
gested as a clue towards a resolution of the 'missing satel- 
lite problem' of cold DM simulations [21] . In [22] it was 
found that once a sub-halo falls into a larger halo mass- 
loss occurs continually through tidal stripping and the or- 
bit of the sub-halo decays down towards the centre of the 
larger halo. Further it was found that the mass-loss can 
be significant, resulting in only ~ 2% (8%) of the mass 
of a sub-halo accreted at z = 2 (1) surviving, and this re- 
sult appears to be independent of the masses of the halo 
and subhalo. The tidal stripping of DM from primor- 
dial globular clusters has been investigated with several 
N-body simulations (see e.g. [23H1S])- Results suggest 
that globular clusters can be formed naturally within DM 
subhalos which are subsequently tidally stripped by the 
host galaxy, resulting in a baryon-dominated core with 
a small mass-to-light ratio, resembling observed globular 
clusters. In particular a recent analysis of the Aquar- 
ius simulation [26j lends support to this scenario, and 
an approximate relation between the current mass of a 
globular cluster and the mass of the initial subhalo it was 
embedded within is given as Mac — 0.0038 MoMfi- 

A recent review [37] also argues that metal-poor glob- 
ular clusters formed in low-mass DM halos in the early 
universe. 

The observed cold WDs reside in the dense core of M4, 
which has survived previous tidal stripping until now. 
Therefore it is reasonable to assume that the majority 
of the DM in the core of M4, well within the tidal ra- 
dius, will also have survived from the early subhalo. This 
assumption is supported by the results of [23' where it 
is found that the presence of the stellar core makes the 
subhalos more resilient to tidal stripping, and for NFW 
subhalo profiles the DM density in the innermost regions 
of the subhalo is not modified by the external tidal field. 
Outside of the star dominated region the DM subhalo 
is stripped back to the tidal radius, thus resulting in a 
mass-to-light ratio close to the purely baryonic value. 

Similar reasoning has led to recent consideration of in- 
direct DM signals from DM annihilation in other globu- 
lar clusters, [2B] and the VERITAS collaboration argue 
that the association of globular clusters and DM halos 
fits naturally into the standard paradigm of hierarchical 
structure formation [5S]. 

Motivated by these developments we follow similar 
methods to those used in [21, to which we refer the 



3 



reader for details. The mass of baryonic matter in M4 
is estimated to be ^ lO^M© and the core radius of 
0.83' in arcminutes imphes — 0.531 pc when com- 
bined with a distance to the cluster of 2.2 kpc. The tidal 
radius is estimated using a concentration parameter of 
log(rt/rc) — 1.59 giving rt — 20.66 pc. These parame- 
ters set the baryon density distribution, which we model 
with a King profile. 

Using cosmological data and taking mass loss during 
stellar evolution into account, the amount of DM in the 
original M4 subhalo is estimated to be Mdm ~ 10 ''M0. 
For details of this estimation see [2 . The virial radius, 
which sets the initial DM distribution is estimated using 
the fitted form of the spherical collapse overdensity [50] : 

^ ^ IStt^ + 82(a„(z) - 1) - 39(a„(z) - 1)^ 

^rn (^) 

where the matter density is given by [31]: 



1 



r!„(i + z)3 



(2) 



We take ilm(O) = 0.273 giving A = 357. The concentra- 
tion of low mass halos is given in [3T] as; 



c(z) 



27 / MpM 

l + z \ WMq 



-0.08 



(3) 



and combining this expression with those for the virial ra- 
dius, scale radius and central density from [2] the original 
DM subhalos are completely determined by the parame- 
ters: 
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0.37 



We model the original DM halo with an NFW profile 
[32]. As discussed in [2 the core density is a very weak 
function of the total mass of the subhalo, changing only 
by a factor 3 for halo masses between IQ^Mq and 10* M0. 

It remains to consider the effects of the baryonic core 
on the DM distribution. Although the DM density may 
be enhanced in the core due to the presence of the bary- 
onic core [33H35] the heating of DM particles due to in- 
teractions with stars may tend to wipe out this enhance- 
ment. Therefore by estimating the timescale over which 
this process occurs with Eqn. 3a. of [3S] we can find the 
radius at which this timescale is equal to the age of the 
universe. We find that this radius lies at r^eat = 1.4 pc 
and, as this is smaller than the radius where the WDs 
are observed, we expect heating effects to be small here. 
The possible important effect is therefore the contraction 
of the DM core due to conservation of angular momen- 
tum when the gas in the original halo which eventually 
forms the globular cluster loses energy and falls into the 
core. We use the algorithm of Gnedin [34j to perform 
this baryonic contraction. Finally as mentioned earlier, 
to take account of the likely tidal stripping of the stars 




0.5 1.0 5.0 10.0 

Radius [pc] 



50.0 



FIG. 1: Densities of stars (solid line) and the DM lialo with 
and without baryonic contraction effects in dashed black and 
dot-dashed red respectively. The region from which the ob- 
served WD data is taken is indicated in shaded blue. The 
vertical dotted lines denote the radius at which heating ef- 
fects start to become important (dotted black) and the tidal 
radius (dotted red). 



and DM halo we truncate the density distribution at the 
tidal radius. 

The estimated DM and star densities for the halos are 
plotted in Fig. [l] One can see that this estimate for the 
DM density is much smaller than the baryonic density, 
and in fact for the contracted (uncontracted) halo DM 
makes up less than 43% (41%) of the total mass of the 
cluster, consistent with the observed lack of DM in glob- 
ular clusters. Further, the total estimated DM content is 
7.7 X IO^Mq, less than 1% of the original 10^ Af© halo. 

We assume a DM density at the largest radius within 
which the WD data is observed, r^ax = 2.3 pc, giving 
Pjjj^.^ = 21 M0 pc"3 = 798 GeV cm"^ for the contracted 
halo [Sg. 



IV. CAPTURE OF IDM IN WHITE DWARVES 

The capture of DM by scattering in stars or planets 
has been studied for some time, see e.g. [^ and 
recently the capture of iDM in the Sun has been studied 
[51-fTU]. It is this work which we extend to include cap- 
ture in WDs and we follow the formalism first set out 
in [8j , which was subsequently extended to include spin- 
dependent scattering as well as spin-independent scatter- 
ing in [XOJ • 

Recently spin-dependent inelastic scattering has been 
suggested as a viable alternative to the standard spin- 
independent iDM scenario [45) . where it is shown that 
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spin-dependent couplings to protons, and not neutrons, 
can give a good fit to the DAMA data whilst remain- 
ing consistent with other experiments. This scenario is 
subject to limits from scattering in the Sun [TU] however 
limits from capture in WDs should be weak as WDs are 
mostly composed of ^^C and ^^O which have no nuclear 
spin. It is tempting to consider limits from scattering on 
^■^C which makes up ~ 1% of all carbon and has nuclear 
spin / — 1/2, however this spin is carried by an unpaired 
neutron and thus the scenario described in [45] will lead 
to negligible capture rates in WDs. Therefore we only 
consider spin-independent scattering iDM in this work. 

To calculate the capture rate in a WD we use the equa- 
tions contained in Section II. of [TU]. To take account of 
the incoherent scattering of DM in the nucleus we also 
use the Helm form factor [56] and follow [2] in making 
the conservative assumption that the WDs are entirely 
composed of carbon. We also assume the DM couplings 
to neutrons and protons are the same. This can be cor- 
rected for specific models by re-scaling the cross-section 
accordingly. 

To find the DM velocity dispersion we make the as- 
sumption of hydrostatic equilibrium and integrate the 
hydrostatic equation for spherical geometry using the 
baryon and DM distributions shown in Fig. [Tj We find 
that the DM velocity dispersion doesn't exceed 8 km s~^ 
and, as the capture rate decreases with increasing disper- 
sion, we set Vo to this value. Similarly the WD velocity 
through the DM is likely to be of the order of the veloc- 
ity dispersion and we find this doesn't exceed 6 km s~^, 
however to set conservative limits we set v^, = 20 km s~^ 
which is the escape velocity at the inner radius at which 
the WDs are observed. 

We calculate the escape velocity and density of nuclei 
within a given WD using the Salpeter equation of state 
[T7] . Due to the large escape velocity of a WD the typ- 
ical kinetic energy of an in-falling DM particle is of the 
order ~ 1 MeV. Therefore all of the in-falling particles 
easily have enough kinetic energy to overcome the inelas- 
tic splitting and scatter. This makes the inelastic split- 
ting relatively unimportant up to splittings S 1 MeV, 
where the capture rate starts to decrease. This is shown 
in Fig. [2] Although the splittings associated with iDM 
are much too small to decrease the capture rate signifi- 
cantly we include them in our calculations for the sake 
of thoroughness. 

The capture rate typically falls as the inverse of the 
DM mass, therefore the luminosity should be largely in- 
dependent of the DM mass. However there is a subtle in- 
terplay between two factors which leads to a dependence 
not only on the DM mass but also on the WD mass. 
The first factor is due to the conversion between a DM- 
nucleon and DM-nucleus cross-section at zero momentum 
transfer, which results in a factor of a^pf oc {fi^pf / fi^n)'^ 
where fi is the reduced mass and N (n) subscripts denote 
the nucleus (nucleon). This factor has a preference for 
heavy DM particles. However there is also suppression 
due to the nuclear form factor [371. Therefore, although 
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FIG. 2: The DM capture rate for a solar mass WD as a 
function of inelastic splitting 5. The DM parameters are = 



50 GeV and (t„ 
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om . The capture rate is largely 



independent of the inelastic splitting up to J ~ 1 MeV, when 
it starts to fall off rapidly. 



heavier DM particles have a larger range of scattering 
energies, the higher energy events are suppressed. This 
effect therefore discriminates against heavy DM particles. 

Which of these two factors wins out depends on the 
WD mass. As heavy WDs have greater escape veloci- 
ties (~ (7 - 12) X 10^ km s~^), heavy DM particles feel 
the form factor suppression more and light DM particles 
lead to greater luminosities. For light WDs the escape 
velocities are lower (~ (2 — 3) x lO'^ km s~^) the form 
factor suppression is subdominant and heavier DM par- 
ticles lead to a greater luminosity. The mass dependence 
for two different WDs is illustrated in Fig. |3] 

As described in [5] the time-scale for equilibrium be- 
tween the capture and annihilation of DM in WDs is 
roughly of the order of one year, and therefore we can 
safely assume that the capture rate is one- half of the anni- 
hilation rate. We also assume that all of the energy of the 
annihilating particles contributes to the black-body lumi- 
nosity of the WDs, however specific DM models could be 
investigated by calculating the fraction of energy lost as 
neutrinos per annihilation and weakening the limits ac- 
cordingly. 



V. UNCERTAINTIES 

We will now discuss some of the assumptions that have 
gone into this calculation. The greatest source of uncer- 
tainty is the estimation of the DM density. It should be 
noted that there are some models where globular clusters 
are not formed due to the collapse of a DM dominated 
halo, in particular the observations made by Gilmore and 
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FIG. 3: The luminosity of a WD due to capture and anni- 
hilation of DM as a function of DM mass for WDs of mass 
= Mq (solid line) and M* = 0.2 Mq (dashed). The DM 
parameters are 5 = 130 keV and a„ — 10~*^ cm^. There is a 
slight enhancement around ~ 10 GeV for the solar mass 
WD and the luminosity is suppressed for small DM masses. 
This behavior is discussed in the text. 



collaborators have lead them to argue that low mass star 
clusters are fundamentally different to higher mass galax- 
ies rather than both being members of a continuous fam- 
ily [46'. The explanation for this scenario typically re- 
quires some kind of modification of DM such as warm 
DM, a cold/hot admixture or a non-zero self-interaction 
cross section such that there is a minimum size for dark 
matter halos in the Universe. Since we are trying to put 
constraints on models of cold DM, it is a consistent as- 
sumption that the globular clusters do form in the centre 
of DM halos in the early Universe but we note that this 
is an uncertainty. 

It has been shown that for direct detection experiments 
the details of the DM velocity distribution can have a 
significant impact on detection rates [47] , particularly for 
iDM [JSJHS]. However, due to the large escape velocity 
of the WDs all in-falling DM particles will have a large 
enough velocity to scatter and the details of the velocity 
distribution will be unimportant. 

We have made the estimate that the WDs are travel- 
ing at the local escape velocity, however a more realistic 
(but less conservative) assumption would be that they 
are traveling at a speed closer to the local velocity dis- 
persion, which is roughly a factor of 3 smaller. As the 
capture rate is inversely proportional to this speed we 
may have underestimated the capture rate by the same 
amount. 

The observed WDs may not be entirely composed of 
Carbon, and may contain heavier elements, however this 
assumption is safe as the capture rate is smallest for light 
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FIG. 4; Observed WDs (black dots) and main-sequance stars 
(open red circles). Also plotted are the luminosity and tem- 
perature of WDs in the mass range 0.1 — 1.35 Mq for the two 
benchmark points il and 12 (dashed) described in the text. 
For both curves a DM-nucleon cross-section of 10~^^ cm'^ is 
assumed. The luminosity from DM capture is greater than a 
significant number of the observed WD luminosities. 

target nuclei. Due to the large escape velocity the energy 
transfer in scattering events can be large (~ 1 MeV), 
and the scattering is thus significantly suppressed by the 
nuclear form factor. We use the Helm form factor, which 
for light nuclei can be up to ~ 20% greater than more 
realistic form factors at these high energy transfers |50j . 
Therefore a conservative estimate of the uncertainty due 
to the choice of form factor is of the order 20%. 

For considerations relating to errors in WD observa- 
tions we refer the reader to |15) . 

VI. RESULTS 

In Fig. |4] we show the observed WDs in the 
Temperature-Luminosity plane. On the same plot we 
show curves for WDs whose sole energy source is due to 
DM annihilation in the core for WDs ranging in mass 
from 0.1 Mq to 1.35 Mq. These curves correspond to 
two benchmark points: 
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The mass and splitting for the point il corresponds 
roughly to best-fit points for channeled scattering on Io- 
dine and in [45l El] it is shown that for these parameters 
consistency with other direct detection experiments is 
possible. The mass and splitting in point 12 corresponds 



6 



to the conventional quenched scattering iDM scenario. 
This scenario can be considered ruled out by other direct 
detection experiments, however this depends sensitively 
on how (or whether) data from different experimental 
runs is combined, the inclusion of borderline scattering 
events in the CRESST detector, choice of Iodine quench- 
ing factor and details of the DM velocity distribution 
[48. 49, 51 . Therefore for completeness we still include 
this choice of parameters in our analysis. We choose a 
cross-section of an = 10"^"'^ cm^ as this is below the cross- 
section at which the optically thick limit applies and the 
capture rate becomes independent of the scattering cross- 
section. This important point was first emphasized in the 
context of inelastic dark matter capture in [52] which ap- 
peared shortly after the initial preprint of the current 
work was placed on the ArXiv. We have updated our 
cross-sections so as to stay below this limit although the 
main conclusions of this paper are unchanged. 

As the temperature of a black-body is related to the 
luminosity as T oc L^/^ Fig. [4] can be misleading, as a 
change in cross-section does not correspond to a simple 
re-scaling of the luminosity. Therefore in Fig. [5] we esti- 
mate the observed WD masses as described in Section. 
[TTl and plot curves showing the luminosity due to DM 
capture and annihilation for a given WD mass. 
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FIG. 5: Observed WDs (black dots). Also plotted are the 
luminosity and temperature of WDs in the mass range 0.1 — 
1.35 Mq for the two benchmark points il and i2 (dashed). 
The upper red curves correspond to a cross-section of cr„ = 
10~^^ cm^ and the lower black curves to (t„ = 10"'*^ cm^. 

As can be seen the low luminosities and temperatures 
of the observed WDs appear incompatible with cross sec- 
tions greater than ct„ ~ 10"''^ cm^ for either of the 
benchmark points. A cross-section below 10^^^ cm^ 
could possibly be argued as acceptable given the uncer- 
tainties, and a cross-section as low as lO""*^ cm^ would 
evade these bounds entirely. However any cross-section 



greater than 10~^^ cm^ would certainly appear to con- 
tradict the existence of such cold WDs in M4. 

Recent analyses of iDM find a best fit cross-section of 
2 X 10-37 cm^ [ig and 4.7 x IQ-^^ cm^ [sT for the bench- 
mark point il. As these cross-sections are greater than 
^ 10~*^ cm^ both scenarios appear excluded by the anal- 
ysis above, and even the lowest cross-section found for the 
channeled iodine region of 9 x 10^''^ cm^ (51j would be 
difficult to reconcile with these results. The conventional 
iDM scenario of quenched scattering off iodine (bench- 
mark point 12) typically requires cross-sections greater 
than ~ 3 X lO"'"' cm^, with a best-fit point at cr„ = 10"^* 
cm^ [5S]. This scenario seems excluded by this analysis. 

These constraints can be evaded in any combination of 
the following scenarios: 

• The DM density is less than ^ 1% of that estimated 
here. This appears plausible, however it would im- 
ply that DM makes up less than ^ 1% of the total 
mass of the globular cluster. If globular clusters 
are born in dark matter halos it is hard to imag- 
ine how the ratio of DM to baryonic matter in M4 
could be so far below cosmological values. If M4 did 
not form in a DM halo it is likely the DM density 
would be low enough to evade these bounds. 

• The DM could annihilate predominantly to neutri- 
nos, thus contributing little to the visible luminos- 
ity of the WDs. This scenario would be difficult to 
reconcile with limits on neutrino fluxes from DM 
annihilation in the Sun [8UIO). 

• The iDM couples to nuclei only through spin- 
dependent interactions as recently suggested |45) . 
If iDM couples to neutrons then it may be possible 
to set limits by considering scattering off ^"^C, how- 
ever as this scenario already appears disfavoured 
by direct detection experiments |45j this is not in- 
vestigated here. If the coupling is only to protons 
then limits from WDs pose little threat. 

From simple assumptions about the DM density in 
globular clusters (including their formation in DM ha- 
los), the composition of cold WDs, and the capture of 
iDM we argue that the iDM explanation of the annual 
modulation observed by the DAMA collaboration is in- 
compatible with the observation of cold WDs in M4 if 
the DM annihilates to standard model particles. Alter- 
natively if the inelastic dark matter scenario were to be 
confirmed and it was found to annihilate to standard- 
model particles then this would imply a much lower dark 
matter density in the core of M4 than would be expected 
if it were to have formed in a dark matter halo. 



VII. DISCUSSION 

We now discuss some of the salient features of white 
dwarves which make them a unique probe of DM. 



7 



It is interesting to note that in the case of elastic DM 
scattering cross-sections (t„ ~ 10""*^ cm^ evade the WD 
bounds [S5]. It is unhkely that observations of WDs much 
cooler than those in M4 will be made as the low lumi- 
nosity cut-off has been observed, and the luminosity of 
WDs is limited by the age of the Universe. Therefore it 
is unlikely that limits from WDs will ever compete with 
direct detection limits for weak-scale elastic DM. 

However WDs constitute unique DM probes for three 
reasons: 

• The large escape velocity enables in-falling DM 
particles to easily overcome inelastic splittings and 
leads to large energy transfers in scattering. 

• The low mass of carbon gives WDs sensitivity to 
light DM scenarios, where most direct detection ex- 
periments lose sensitivity. 

• Limits from capture in the Sun arise due to neu- 
trino annihilation products and are therefore insen- 
sitive to DM annihilating to e~^e~ , fi'^fJ,~, 77, light 
hadrons or gluons. It is specifically this scenario 
where limits from WDs are strongest. 

We have only considered iDM capture in this work. 



however numerous possibilities exist for future study of 
DM capture in WDs. Examples would include DM with 
mass splittings of the order a few MeV or DM which scat- 
ters through a light mediator, ^ MeV, which could be 
enhanced in WDs through the propagator — 

Finally we note that if DM were to be discovered in 
future experiments, and details of the DM-nucleon cross- 
section and annihilation products were to be established, 
then cold WDs could be used to determine an upper limit 
on the DM density within M4, thus giving clues as to the 
formation of globular clusters. 
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